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The interaction between mechanical devices and light has been the focus of a very 
intensive research activity, in the last six years, directed towards the demonstration of 
quantum control achieved and operated over quasi-macroscopic systems undergoing quite 
non-trivial dynamics [IJ. This has culminated in the very recent experimental demonstrations 
of ground-state cooling of massive mechanical modes achieved through radiation pressure- 
induced dynamical back-action [2J, the realization of important technological steps directed 
towards the construction of quantum-limited position sensors and transducers [3], as well 
as the engineering of non-linearities strong enough to mimic, in a fully mechanical 
context, processes and dynamics that were so far typical of quantum optical settings BUl. 
The perspective of exploiting massive mechanical systems in order to prepare quantum 
states and enforce non-classical dynamics has spurred a few theoretical proposals for 
the achievement of entanglement involving both optical and mechanical systems (from 
here on referred to as optomechanical entanglement for brevity) [j5J. Such a scenario 
has been enlarged to consider schemes for the generation of fully mechanical entangled 
states [[6l |7l| and quantum correlations in hybrid systems involving light, mechanical 
modes and collective atomic ones [8]. Such achievements are endowed with intrinsic 
fundamental interest as they open up the way to the experimental study of mechanical 
versions of the long-craved Schrddinger cat and Einstein-Podolski-Rosen state, as well as 
the investigation of exotic decoherence mechanisms including the quantum gravitational 
collapse model [10], the continuous spontaneous localization process [1 IJ and gravitationally 
induced decoherence fil2J . Moreover, the considerable flexibility and high potential for 
hybridization of optomechanical systems leave suflicient space for a significant re-design 
of the way a device for quantum communication is thought of. Taking inspiration from 
the (by now well grounded) paradigm of quantum networks for communication made out 
of local processing nodes, mutually interconnected by quantum channels [fT3l [T4l . one can 
think about a web of remote and non-interacting optomechanical systems that share multi- 
photon light fields ensuring long-haul connectivity. In a sense, the embryo of such an 
architecture was already in the work by Zhang et al. [ITSl . where the coupling between a 
two-mode squeezed vacuum (TMSV) state and two independent optomechanical cavities was 
explored for mechanical Hamiltonian engineering. This scheme has been pushed further 
in BT], where an extensive study of the efficiency of entanglement distribution from purely 
optical TMSV resources to all-mechanical receivers, including the proposal of a mechanism 
for the measurement of such entanglement, was presented. It should be noticed that alternative 
configurations for the distribution of entanglement in a mechanical network itself have been 
proposed, based on optical postselection [|T6l or the use of atomic information carriers 
propagating in close proximity of monolithic microresonators [ITtI . 

In this paper, we go significantly beyond the proof-of-principle study reported in [|71 
and analyze the performance of mechanical entanglement distribution by means of more 
general multimode optical resources to determine the extents of the validity of such 
paradigm. We start by digging deep into the two-mode resource case, which entails the 
minimal configuration for entanglement distribution and, as such, the basic building block 
for any optomechanical network. By working with generally mixed two-mode Gaussian 
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states, we show that the mechanism is very much relying on the symmetry of the optical 
resource's state: small asymmetries in the statistics of the two modes composing the optical 
entanglement carrier correspond to a spoiled mechanical entanglement distribution efficiency. 
The distribution mechanism is optimized by the use of TMSV states whose degree of 
squeezing (and hence entanglement) is determined by the working point of the mechanical 
receivers. A posteriori, this demonstrates that the study in fTl entails indeed a configuration 
of optimal entanglement distribution. We then move to the simplest multipartite setting, 
consisting of a (generally mixed) three-mode driving field guiding the dynamics of three 
remote and independent optomechanical cavities. We show that, under conditions of fully 
symmetric multi-mode driving fields, tripartite inseparability is indeed in order as signaled 
by the violation of the positivity of partial transposition criterion [1181 . Our work embodies 
the first quantitative investigation on the influence that the parameters characterizing a multi- 
mode driving field have on the performance of mechanical entanglement distribution. Our 
results allows for the identification of the optimal resources to use at a given optomechanical 
working point as well as, conversely, the way an optomechanical network should be adjusted 
in order to fully exploit the entanglement brought about by the driving field. 

The remainder of this manuscript is organized as follows. In Section [H we introduce 
the model used for entanglement distribution among mechanical receivers that are part of 
simple optomechanical devices and sketch the analytical procedure to determine the stationary 
reduced mechanical state resulting from the interaction with the driving field. Section |2] 
discusses the two-mode case and the corresponding efficiency of the entanglement distribution 
process. This analysis is then extended to a three-site case in Section |3l where symmetric 
and genuinely tripartite entangled fields are used in order to get a three-mode inseparable 
mechanical state. Finally, in Section |4] we draw our conclusions and discuss a few open point 
that we aim at addressing in future. 

As a final remark of this Introduction, it is important to stress that the working points 
used for our simulations adhere perfectly to the current experimental state of the art, while 
the sort of optical states used in order to drive our small optomechanical networks are all well 
within the reach of linear optics experiments aiming at the engineering of continuous-variable 
states. 

1. Steady-state mechanical entanglement from optomechanical dynamics 

We now provide a general description of the mechanical entanglement distribution scheme 
without limiting the number of sites of the network that we consider. Figure [1] shows a pictorial 
representation of the experimental configuration at hand, where the optical circuitry needed 
in order to generate an optical A^-mode entangled resource is condensed in an unspecified 
"Resource generator". On the other hand, the mechanical network consists of remote 
and mutually independent optomechanical resonators. In order to provide a simple explicit 
example, we consider the case of Fabry-Perot cavities each endowed with a light mechanical 
end-mirror, which is one of the configurations currently utilized experimentally [IJ and 
has been exploited to provide the first experimental evidence of radiation pressure-based 
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Figure 1. Mechanical entanglement distribution scheme. A multi-mode state generator 
prepares the continuous variable (CV) resource, which is distributed to the nodes of the 
optomechanical network. The latter consists of remote and non-interacting cavities, each 
endowed with a light movable mirror of mass effective m, mechanical quality factor Q and 
natural frequency w,„. Each resonator leaks photons at a rate k, has length L, free frequency 
(jjc and is driven by a quasi-resonant classical pump. The movable mirrors are at equilibrium 
with their respective environments at temperature T . 



passive cooling [fT9l and demonstrate strong optomechanical coupling [|20l . Throughout the 
manuscript, for the sake of simplicity, we assume identical optomechanical nodes. Therefore, 
each cavity has length L, natural frequency ojc, photon decay rate k and is terminated 
by a mechanically movable mirror of effective mass m, free frequency (under harmonic 
approximations) cu,,, and mechanical quality factor Q. Cavity j = 1,2, .., A'^ is driven by both 
a large-intensity classical field (pump power rj and frequency coo) and mode Mj of the (weak) 
optical entangled resource. The drivings populate the cavity field, which in turn is coupled, 
via radiation pressure at bare rate x = <^c/^ [21] . to the mechanical mirror. The free energy 
of the whole optomechanical network thus reads 

N N 



-xqj)a]aj + -^p] + ]^ma?„q] + TiEyj 



(1) 



where 6 = a>, -ojo is the cavity-pump detuning, E = y ^ is the amplitude of the displacement 
undergone by the cavity field due to the strong pump, a (a'^) is the corresponding 
annihilation (creation) operator [with yj = i(a J - aj) / V2 the out-of-phase quadrature of 
the cavity field] and {qj^Pj} the position and momentum operator of the mechanical 
oscillator satisfying the canonical commutation relations [qj,Pk] = ihSjk ij,k = l,..,N). 
The dynamics induced by this Hamiltonian is captured by the quantum Langevin equations 
of the whole optomechanical system (we introduce the vector of operators for subsystem j 
Oj = (qj,Pj, Xj,yj) with xj = {a] + aj)/ V2) 

dtOi = ^[H, Or] + noise terms ( / = 1, .., A^), (2) 
n 

which incorporate the mechanical damping at rate 7,,, = a),„/Q, the cavity photon decay and 
two noise contributions that deserve a more detailed introduction. The first is the incoherent 
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Brownian dynamics induced by the interaction between the mechanical mode and the 
ensemble of phononic modes excited in the substrate onto which the mechanical oscillator 
is fabricated [|22]. In principle, this mechanism gives rise to non-Markovian evolutions, as 
encompassed by the two-time correlation function of the (zero-mean) Langevin force operator 
bj describing Brownian motion [|23l 



{bj{t)bj(t')) = J 



, , fico 

coth + 1 

2k^T 



dco. (3) 



Here T is the temperature of the phononic background modes (with Ohmic spectral density) 
with which each mechanical oscillator is assumed to be in thermal equilibrium and Icq is 
the Boltzmann constant. Evidently, the noise affecting the mechanical mode is, in general, 
not delta-correlated. However, for large mechanical quality factors, {bj{t)bj{t')) becomes 
approximately delta-correlated, thus reinstating a Markovian description of the associated 
noise. The validity of such approximation is guaranteed across this paper. 

The second source of noise is due to the open nature of the cavities, which determines the 
effect of the weak entangled driving fields onto the dynamics of the mechanical oscillators. In 
fact, in our model, the pumped optomechanical network is affected by an A'^-mode quantum- 
correlated environment that we use in order to distribute mechanical entanglement. The details 
of the correlations shared by the environmental modes depend on the form of the state that we 
decide to use as an entangled resource for the distribution mechanism and will be specified 
later on in this paper by assigning the appropriate vector of two-time correlation functions 
Cin{t, t') (with t and t' two arbitrary instants of time in the system's evolution). 

Finally, under the assumption of intense pump fields, we can focus on the dynamics of 
the quantum fluctuations 60 j around the mean value of each operator Oj, which are thus 
expanded as 

dj^0j + 6dj (j=\,..,N). (4) 

Each Oj that is determined by the steady-state Langevin equations [|22]| . It should be stressed 
that the optomechanical subsystems are mutually independent, so that each light-mechanical 
mode interaction occurs regardless of similar events taking place at other sites of the network. 
In this sense, by neglecting the sources of noise, the time-evolution operator of the whole 
system would read U(t) = ®^^^e~^"j' with Hj that rules the coupling between mode Mj and 
oscillator j, as specified in ([T])- The quantum correlations shared by the input noise link the 
various optomechanical subsystems to each other. 

While details on specific and interesting instances of such correlated environments will 
be given in the next Sections, here we briefly sketch the way one can track down the evolution 
of the system. First, the Langevin equations for the quantum fluctuations introduced above 
can be compactly cast into the vector form dt[6dj(t)] = KjdOjit) + Nj{t) with the dynamical 
matrix 

1/m 

-mo^l -r„, V2%2l(c,J V2%3(c,,) 
-V2;t-3(Css) A 
V2;t2l(Css) -A -K 



Kj = 



(5) 
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and the vector of noise fluctuation operators Nj = (0, bj(t), yJlK5x^J{t), y[2K6fJ{t)), which 
depends on the Brownian-motion Langevin force and the quadrature fluctuation operators of 
the input mode j = 1, .., A^. We have introduced the steady-state amplitude of the field in 
each cavity Css = EI{k + iA) and the detuning A = 5 - fix^\CiJ^ K^oj^m)- It should be noted 
that, notwithstanding the non-linear character of the light-mechanical mode interaction, the 
description in terms of quantum fluctuation operators gives back a fully linear model that 
guarantees the solvability of the problem encompassed by (|5]) ll5l 12211711. 

Here we will be interested in the steady-state regime achieved at long interaction times 
after all transients induced by the light-mechanical mode interaction have faded away. We thus 
drop the time derivative in the Langevin equations and solve the associated time-dependent 
set of linear algebraic equations in the Fourier domain. In the following, we will concentrate 
on the properties of the mechanical system only. We thus report the solution of the Langevin 
equations for the position and momentum operators of the mechanical oscillator, which 
read 

moj{ ^|TKflxi.Cs.5C'f{-^ + iK + oS) + 5C'''c*(^ + i/c + oj>)) + i5Bi^^ + {k- ioof)] 

S'P ■ = ^ — — (6^ 

' 2|c,,PMa'2 + miiy^co + co^ - coi) (A^ + (k - ioj)^) 

i ^^2i^nx[c,s^cY'(-^ + iK + co) + 6Cfc*JA + iK + oj)] - 6Sj[A^ + {k- itof] 

In the equations above, 50j stands for the Fourier-transformed fluctuation operator 60 j (the 
dependence on frequency is understood) and 6k'f = {6c'f+dcJ)/ V2, 6yf = i{5c'f -di'f)! V2. 
Quite clearly, the dynamics of the mechanical oscillator depends crucially on the input 
noise fluctuation operators: from © and (|7]) it is evident that any correlation function of 
the mechanical variables will be determined by the input noise correlators and, in turn, the 
entanglement that the engineered environment brings about. In the next Sections we specify 
the set of input correlators and determine the entanglement that is distributed among the 
elements of the small mechanical network that we are working on. Finally, throughout our 
simulations we will use the set of parameters for the optomechanical system given in Table HI 
which match the current experimental state of the art, thus making our proposal very close to 



Table 1. Parameters used for the simulations run throughout the manuscript (we assume the 
same values for each optomechanical subsystem). The set of parameters is taken from the 
experiment in Ref. [201, where a slightly lower input power was used. 



Parameter 


Symbol 


Value 


Mechanical mass 


m 


145 ng 


Mechanical frequency 


tojln 


947 KHz 


Cavity length 


L 


25 mm 


Input power 




20 mW 


Cavity-field wavelength 


A 


1064 nm 


Optical damping rate 


Kjln 


215 KHz 


Mechanical quality factor 


Q 


7000 
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feasibility. 



2. Two-site optomechanical network 

Here we restrict our study to the case of N = 2 sites of the network, which is thus driven by 
a two-mode quantum correlated state of light. The linearity of the dynamics of the quantum 
fluctuations ensures that the corresponding map preserves the Gaussian nature of an input 
optomechanical state. This observation allows us to limit our attention to just this class 
of continuous-variable states. In fact, Gaussian states are those most easily generated in 
quantum optics laboratories. Moreover, Gaussian resources appear to maximize the efficiency 
of schemes for the transfer of entanglement from continuous to discrete-variable systems, 
when compared to experimentally accessible non-Gaussian states such as photon- subtracted 
or added fields [14J. Finally, as a key assumption of our model is that each mechanical 
oscillator is initially in thermal equilibrium with its own background of phononic modes, 
we take the elements of the mechanical network as initialized in thermal states with mean 
occupation number nj = [e*B'^ - 1]" . The Gaussian-preserving nature of the problem allows 
us to use the formalism of covariance matrices: given an A'^-mode systems described by the 
2A'^-long set of quadrature operators Q = {qi,pu ..,qN, Pn), we define the covariance matrix 
elements (cr)^; = (QtQi + QiQk)/^ {k,l = 1,2, ..,2N). The covariance matrix of a physical 
state should satisfy the Heisenberg-Robertson uncertainty principle cr -i- iQ,f^/2 > with 
= ©^licTy the so-called symplectic matrix and o-y the y-Pauli matrix. The Gaussian state 
of the system is determined by the assignment of cr. 

We now need a general description of the two-mode resource that we are going to employ. 
The covariance matrix of any (pure or mixed) two-mode Gaussian state can be put in the so- 
called standard form parameterized as 



1 

"=2 



a c+ 

a c 

c+ b 

c- b 



(8) 



with 



a = s + d, b = s - d, c+ = _ (9) 

4y/s^-d^ 

and /v = 4x^ + (g^ + l)(A - l)/2 (x = d, s). The state corresponding to cr is physical and 
entangled for s > l,\d\ < s - l,2\d\ + I < g < 2s - I and A e [-1, 1]. Here, s and d set 
the local entropy of each mode, g determines the global purity of the state, while the degree 
of entanglement shared by the two modes grows with A. The assignment of the input optical 
covariance matrix cTopt specifies the set of two-times correlators Cin(t, t') needed to solve the 
mechanical dynamics. Indeed, a covariance matrix such as ^ is obtained by taking an optical 
state having the only non-zero correlators as 

Qnit, t') = i{6cf(t)6t^it')), <5c;"(05cf (?')), {6cf(t)6c'^(t')}, (10) 
<54"(05£f (?')), {Sc';'(t)6c^(t')), {6c'f(t)6c'^Ht'))) (11) 
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Jci>,„(t-t') 



,,)(C++C-) 



6{t - 1'). 



(12) 



2 2 

Physically, this corresponds to taking a weak driving field whose modes have frequency 
0*0 + Oi>m [7]. The covariance matrix cTj^ech resulting from the assignment of a cTopj equal 
to ([8]) maintains the standard form and reads 



""mech 



{ m]) 





{SQ16Q2) 








{5P]) 





{SP16P2) 







(SQl) 





\ 


{6P16P2) 





(SPj) 



(13) 



where each element is calculated using (|6l) and and then transforming back to the 
time domain [7|. We have generated a random sample of input two-mode states by 
taking parameters uniformly distributed across their ranges of variation, which are in turn 
determined by setting an upper bound to the possible values of s. Correspondingly, we have 
calculated the output entanglement distributed between the mechanical oscillators. This has 
been quantified using logarithmic negativity [|25l . which for Gaussian states is defined as 
Xat = max[0, -ln2v_], where v_ = mineigliQao^l is the minimum symplectic eigenvalue 
of the covariance matrix cr = Pi\2crPi\2 associated with the partially transposed state of the 
mechanical network [24] (Pip = diag[l, 1, 1,-1] is the matrix that, by inverting the sign 
of momentum of the mechanical oscillator 2, realizes partial transposition at the level of 
covariance matrices [26J). The results are displayed in Figure |2] (a), where we compare the 
input optical entanglement (determined calculating the logarithmic negativity of the two-mode 
input state) and the corresponding output mechanical entanglement associated with (fT3]) and 
calculated at the detuning A = co,„ that maximizes it. Although the range of entanglement 
of the input states extends uniformly from to more than 3, mechanical entanglement can 
only be distributed when the input one is larger than a non-zero threshold determined by the 
working parameters at the optomechanical level and not the input field's properties. We have 
also found that the system favors symmetric input states: the number of entangled mechanical 
states obtained when considering input states characterized by J = is much larger than the 
analogous number for asymmetric states. This results from the uniformity of the parameters 
within the optomechanical network and is an effect that can be cured by matching the optical 
asymmetry with an equal mechanical one. Most importantly, the entangled distributed in the 
mechanical network is constrained to stay below a curve that maximizes the degree of output 
entanglement at assigned values of the input one. We have found that the optical resource 
giving rise to such boundary should have d = and g = A = 1, i.e. it should be a pure, 
symmetric state (thus equivalent to a TMSV) with covariance matrix 

V^2_ 1 



Copt - T 



s 








s 




-Vs^ 




s 






s 



1 



(14) 



The behavior of the boundary suggested by the numerical exploration is confirmed by an exact 
analytical approach which we sketch here. First, we write the input optical entanglement e in 
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(a) (b) (c) 




Figure 2. (a) Output mechanical entanglement against input optical one for a symmetric 
two-mode resource with identical single-mode variances (d = 0) and variable degree 
of entanglement. The inner red dots represent the distributed mechanical entanglement 
corresponding to 10'* randomly taken mixed input states. On the other hand, the black points 
on the boundary stand for the mechanical entanglement achieved when 5000 pure = 1) 
symmetric input states (picked up randomly) are used, (b) The numerical bound that is 
highlighted in panel (a) enjoys an analytical form (whose derivation is sketched in the body 
of the paper) which is displayed for comparison with the numerics, (c) We study v_ and its 
derivative with respect to s for a pure, symmetric two-mode resource state. We find v_<l/2 
for values of s up to the the dashed vertical lines. At s = 2.501, v- has a minimum (the 
entanglement has a maximum). 



terms of the parameter s characterizing the optimal resource as 

6 = In 1 1 - + 25 V52 - 1 1 ^ 5 = cosh e. (15) 

This expression for s can then be replaced in the analytical one for v-, the minimum 
symplectic eigenvalue of the partially transposed covariance matrix of the mechanical 
network. While it is possible to get an exact closed expression for v_ as a function of all the 
parameters of the optomechanical setting, this is too lengthy to be reported here. In Figure |2] 
(b) we show such analytic bound for the parameters used in panel (a). We defer an analysis 
of the behavior of such expression against some of the key parameters of our system to later. 

Therefore, a TMSV embodies the ideal resource for the distribution of entanglement in 
the mechanical network, thus proving that the analysis in [jTJ is optimal. Moreover, we can 
determine the value of s that, at the optimal detuning A = co^, maximizes the transferred 
entanglement, thus identifying the optimal input resource for the distribution scheme. To do 
this, we look for the solution of the equation dsV-(s) = at g = A = 1 and for a symmetric 
resource. Figure |2] (c) shows the trend followed by both v_ and its derivative with respect 
to s. At s = 2.501, the symplectic eigenvalue of the partially transposed mechanical state 
achieves its minimum: again, the analytic expression for such solution is too uninformative to 
be shown here. 

While, in line with expectations, we have checked that an increasing mass and/or 
mechanical damping rate reduce the amount of entanglement between the mechanical 
systems, the effects of the deviations from ideality of the resource are yet to be assessed. This 
is relevant given that it is likely that the optical resource produced in a quantum laboratory 
would be mixed (due, for instance, to residual thermal nature at the two-mode state generation 
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(b) 
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Figure 3. (a) Minimum symplectic eigenvalue of the partially transposed covariance matrix of 
the mechanical system studied against the single-mode variance s and the global state purity 
g of the input resource, for d - Q and A - \. Whenever v_ > 1/2, mechanical entanglement 
disappears, (b) The function v-(g, s) is asymmetric: the larger g, the larger the value of s at 
which entanglement is maximized. The range of s within which entanglement is distributed 
to the mechanical network is maximum for ^ = 1 (lower red boundary) and shrinks as g 
grows. Ai g - 5.01 (corresponding to the upper green boundary) all mechanical entanglement 
disappears. 



Stage). We thus now address the effects that a non-unit purity of the two-mode optical resource 
has on the performance of our scheme. In order to do this we have considered the range of 
values of {s, g) for a symmetric optical input state (with A= I) that leaves v_ < 1 /2. As shown 
in Figure |3](a), the scheme tolerates a considerable degree of mixedness of the input resource: 
only at g = 5.01 (corresponding to a state purity of 0.2) and for the optomechanical parameters 
stated in Table [H the distributed entanglement disappears. The value of s that optimizes 
the entanglement distribution increases with g, making the function V-{g, s) asymmetric with 
respect to s. This is seen in Figure |3] (b), where we highlight the region of the {g, s) space 
within which mechanical entanglement is found. 

Another important parameter to consider is the initial temperature T of the mechanical 
systems. How warm should the network be in order to jeopardize the entanglement 
distribution mechanism? The process appears to be robust against thermal effects. Figure|4](a) 
shows that the decrease of the maximum value of the output mechanical entanglement is quasi 
exponential and very flat at small temperatures (the drop in the maximum entanglement from 
T = 10~^K to lO'^^K is only quite marginal. Notice that the horizontal axis is in logarithmic 
scale). Entanglement persists up to T = 0.02K, showing that even macroscopically occupied 
thermal states (at T = 0.0 IK the mean thermal occupation number of a mode at frequency oj^ 
is larger than 10^) can accommodate the inputted entanglement [cf. Figured (a) (inset) and 



(b)]. 



3. Three-site optomechanical network 



We now assume that a three-mode optical driving field is prepared to drive the dynamics 
of a three-site optomechanical network. While the principles of the dynamics at the single 




Figure 4. (a) Output mechanical entanglement against input optical one for input driving 
fields characterized by li = 0, ^ = 1 and A = I. From top to bottom, the curves correspond to 
T = 10 ''K (top-most soUd line), T = 10 '^K (dashed curve), T = lO^^K (dot-dashed curve). 
Inset: we show the behaviour of the maximum output mechanical entanglement against the 
temperature of the mechanical modes (in logarithmic scale), (b) We plot v_ against s for 
^ = 1. The lower (solid) curve corresponds to T = lO^^K, while the upper one shows the case 
corresponding to T = 2.01 x lO^^K. Any situation having T between such extremal values 
gives rise to v_'s falling within gray-shadowed region. 



optomechanical site are the same as those for the two-site case, the description of the optical 
resource needs some care. The covariance matrix of a pure symmetric three-mode state is 
parameterized as [241 



cr = - 



ai 










^13 


^ 





ai 

















a2 





^23 













ai 





''23 


^-13 





^-23 























«3 , 



(16) 



with 



C7i = 



^|\.{ai 


- Ujf - {Uk 




-ajY-{a, + m 






+ Qjf - {flk 


- mm 


+ ajY-{ak + 1)2] 



(17) 
4 ^UiUj 

In order to ensure physicality of the corresponding state, the parameters ay's should satisfy 
the triangular inequality 

- aj\ + l<ak<ai + aj-\ (i j k = \,2,3 with aj > 1). (19) 

The analysis presented in Section |2]has revealed that pure and symmetric optical states enable 
better performances of our scheme, at least for the two-site case. We now conjecture that such 
condition extends to the three-site situation and restrict our attention to the case of aj = a (V j). 
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therefore considering the covariance matrix of the pure three-mode optical state 



'opt 
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x^ 





x^ 
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x~ 
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x+ 
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x+ 








x~ 





a 





x~ 












a 








X' 





x~ 





a 



(20) 



with X- = (a^ - I + Vl - lOa^ + 9a'^^ /(8a). Such a conjecture has been verified a posteriori 
by taking covariance matrices with lower degrees of symmetry (such as cases where the 
variances of only two out of three modes are identical) and checking that the corresponding 
performances of the protocol were inferior to those of the fully symmetric case. The state 
with covariance matrix given by Eq. (|20] ) is genuinely tripartite entangled according to 
the measure given by the so-called minimum residual Gaussian contangle [|27l fires- The 
Gaussian contangle fi^ifo between the pair of subsystems a and b is defined as the convex roof 
(performed all the over pure Gaussian states with covariance matrix c^Qa^^ss strictly smaller that 
the covariance matrix of the state at hand) of the squared logarithmic negativity [|24l and is a 
proper entanglement monotone. For a set of modes fi = {i, j, k, I, .., z}, the Gaussian contangle 
satisfies the monogamy constraint &i\jki..z - Zp={y,*:,/,..,?) Sii// ^ 0^ together with the analogous 
relations obtained by permuting the mode indices at will. With this, for three modes {i, j, k), 
£rps is defined as 



fires = mmn{&i\jk 
n 



fill/t] 



(21) 



with n the operator that permutes the indices [i,i,k}. The residual entanglement of the 
symmetric state with covariance matrix la^^t (the factor 2 is used to uniform the notation 
to the formalism of [|27l ) can be computed analytically as fires = arcsinh^[ Va^ - 1] - 
(I /2)\n^[2a^ + 4-ax~]. Moreover, the state has non null entanglement between any two modes, 
growing (but not indefinitely) with a. Driving fields of this type are within the grasp of current 
experimental possibilities [|24]. 

The corresponding mechanical covariance matrix is the extension of (fT3]) to three 
mechanical modes. Throughout this Section we maintain the very same values of the 
optomechanical parameters given in Table \T\ Needless to say, the output covariance matrix 
Cmech is physical as far as Creech + i^^3/2 > (which is equivalent to the condition on the 
symplectic spectrum mineiglillacri^echl > 1/2). As for entanglement, given the tripartite 
nature of the mechanical network at hand, the situation is more involved than in the two-site 
case. Indeed, we can distinguish the case of bipartite 1-vs-l-mode entanglement, achieved 
by tracing out any one of the three mechanical modes, from the l-vs-2-mode entanglement, 
which we will use to infer the tripartite inseparable nature of the mechanical state. 

Let us start with the first instance, which implies the trace over the degrees of freedom 
of one of the modes of our mechanical network. Needless to say, the enforced symmetry of 
the optical input state guarantees full symmetry of the mechanical one, thus making all the 
two-mode reductions identical to each other. Formally, this is equivalent to consider only the 
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blocks in the mechanical covariance matrix cr^iech containing the degrees of freedom of the 
modes that are not traced out. Entanglement is then calculated by using the very same formal 
apparatus described in Section |2]and we can characterize the output mechanical entanglement 
of the two-mode reduced states versus the analogous quantity for the input optical one [cf. 
Figure |5] (a)] . In analogy with the case of two-site networks, we observe both the non- 
monotonic behavior of the output mechanical entanglement and the existence of a threshold 
in the input optical entanglement for the process to be successful. However, at variance with 
what has been observed in Section |2l the curve of the output entanglement is asymmetric 
with respect to the input one. Needless to say, as the reduced two-site system is mixed, we 
expect the distribution performances in this case to be contained within the optimal boundary 
embodied by a TMSV. This is shown in Figure [5] (b) and its inset. 

Symmetry is fully reinstated when one looks at the l-vs-2-mode splitting of the network. 
Entanglement in such configuration (which is again invariant under permutations of the mode 
indices due to the overall symmetry of the three-mode mechanical covariance matrix) is 
detected by the violation of the Heisenberg-Robertson condition by the partially transposed 
covariance matrix Pi\jk(^mechPi\jk, where Piyk is the matrix that inverts the sign of momentum 
of mode /, thus implementing the partial transposition with respect to such element of 
the mechanical network. Entanglement in such one-vs-two-mode split is then quantified 
as Xa^ = max[0, -ln2v^] with v_ = mineig\iQ.2,P i\jkf^Pi\jk\ (as the partially transposed 
covariance matrix has only one symplectic eigenvalue that can be smaller than 1 /2) and is 
plotted, against the input optical entanglement, in Figure|5](c). Not surprisingly, the similarity 
with Figure l2](b) is only apparent and the one-vs-two-mode entanglement is strictly contained 
within the one-vs-one-mode entanglement. 

With this at hand, we can proceed to the qulitative characterization of the tripartite 
entanglement within the system at hand. Following the classification by Giedke et al. [fTSll . 



(a) (b) (c) 




optical entanglement optical entanglement optical entanglement 



Figure 5. (a) One-vs-one-mode entanglement for a three-site optomeclianical network 
pumped by a pure and symmetric three-mode optical state. We plot the output mechanical 
entanglement against the input optical one for 3000 randomly picked values of a. (b) 
Comparison between one-vs-one-mode entanglement (i.e. the one involving any two modes 
out of the three-site mechanical network) and the two-mode mechanical entanglement as 
determined in Section |2] The first is contained fully within the second, (c) One-vs-two- 
mode entanglement (solid line) for a three-site network for the very same states used in panel 
(a). The curve is fully contained within the area below the optimally distributed two-mode 
mechanical entanglement (dashed line). 
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a tripartite Gaussian state whose one-vs-two-mode bipartitions are inseparable is tripartite 
inseparable. Looking at the results shown in Figure |5](c), this occurs, in our system and for 
the set of optomechanical parameters at hand, for quite a large set of input optical states, thus 
showing that multipartite entanglement used as a resource in our protocol is indeed able to 
generate tripartite entangled states of the mechanical network, in a sense signaling the success 
of the protocol itself. 

Yet, a point that remains to be addressed is the quantification of the genuine tripartite 
entanglement witin the full mechanical state. As the state of the mechanical network is overall 
mixed, the calculation should be performed with care. Without entering unnecessary details 
related to the rather technical approach to the calculation of the residual genuine tripartite 
entanglement in the state of a system of bosonic modes we mention that the results of such 
calculations show that genuine tripartite entanglement is indeed transferred to the mechanical 
network, albeit in an inefficient manner (the average tripartite mechanical entanglement being 
two order of magnitude smaller than the corresponding input one, in general). Yet, this is 
sufficient for the scopes of this investigation, which aimed at providing a formal analysis of 
the distribution performance in the simplest instance of multipartite optomechanical network. 

4. Conclusions 

In this paper we have addressed extensively the problem of entanglement distribution in 
an optomechanical network that is driven by a multipartite entangled optical resource. By 
using the most general parameterization of (generally mixed) Gaussian optical states, we have 
found the conditions under which two-site mechanical entanglement is maximally distributed 
identifying the role that symmetries and local/global state purity of the optical resource play in 
this process. We have demonstrated the efficiency of the protocol, which can be implemented 
using state of the art parameters and the possibility to distribute genuine tripartite mechanical 
entanglement. 

Many questions remain to be addressed, for instance related to the ability of the protocol 
to pass general quantum correlations that go beyond quantum entanglement. This is an 
interesting opportunity in light of the preliminary analysis performed in [7J, where such 
quantum correlations have been shown to be more robust than entanglement to the increase 
of temperature in the mechanical setup. Moreover, it would be interesting to find whether 
optical postselection are able to improve the distribution mechanism, giving rise to more 
substantial genuinely multipartite mechanical entanglement. Finally, we intend to perform a 
classification of the distributed mechanical entanglement at assigned values of the global and 
marginal purity, along the lines of some of the studies in [14], so as to check if a hierarchy of 
mechanical states arises from the distribution process. 
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